Exact Solution for ID Spin-Polarized Fermions with Resonant Interactions 
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Using the asymptotic Bethe Ansatz, we obtain an exact solution of the many-body problem for 
ID spin-polarized fermions with resonant p-wave interactions, taking into account the effects of 
both scattering volume and effective range. Under typical experimental conditions, accounting for 
the effective range, the properties of the system are significantly modified due to the existence of 
"shape" resonances. The excitation spectrum of the considered model has unexpected features, such 
as the inverted position of the particle- and hole-like branches at small momenta, and roton-like 
minima. We find that the frequency of the "breathing" mode in the harmonic trap provides an 
unambiguous signature of the effective range. 

PACS numbers: 03.75.-b, 05.30.Fk 



Experimental progress in the cooling and trapping of 
ultracold atomic gases makes it possible to investigate 
their properties under strong transverse confinement, 
when the motion of atoms is effectively one-dimensional 
(ID). In such experiments the ID interaction parame- 
ters are precisely known and can be tuned using Fesh- 
bach resonances [l[ or by varying the harmonic trans- 
verse confinement strength. Recent experiments @, H| 
have allowed for parameter-free comparison of ID Bosc 
gas properties with a theoretical description based on the 
exactly solvable Lieb-Liniger (LL) model The possi- 
bility to compare experimental results with the outcomes 
of many-body calculations revived an interest in the field 
of exactly solvable ID systems: spin-1/2 fermions [f|, 
Bose-Fermi mixtures 0, 0] , and spinor bosons Q and 
fermions @|. 

In this Letter, we obtain an exact solution for ID 
spin-polarized fermions under resonant scattering condi- 
tions which is relevant for 40 K and 6 Li atoms near 



p-wave Feshbach resonances [111 112} . Such a ID experi- 
mental system has been realized for 40 K 13] . The related 



3D problem has also received significant theoretical at- 
tention recently For spin-polarized fermions only 
scattering in odd partial wave channels is present, and at 
low energies, p-wave scattering is the strongest. If only 
the "scattering volume" is taken into account and the 
"effective range" of p-wave scattering is neglected (see 
Eq. (1) for definitions), then the projection to ID 11511 
results in a fermionic Cheon-Shigehara (CS) model [l6| . 
which is dual to the bosonic LL model. For such a model, 
strongly interacting fermions with resonant interactions 
are mapped to weakly interacting bosons, the so-called 
fermionic Tonks-Girardeau (f TG) limit [17|. However, it 
was shown by L. Pricoupenko [lOj that unlike the case of 
the strongly interacting bosonic TG limit @, 0] , the re- 
quirements for the observation of the fTG limit are quite 
stringent, and the effective range of scattering needs to 



be taken into account. We provide an exact solution that 
accounts for both scattering volume and effective range, 
and obtain significant deviations from the CS model [161 ] 
due to "shape" resonance in the p-wave scattering. We 
find several new effects, such as the inversion of particlc- 
and hole-like spectra for low momenta, roton-like min- 
ima in the excitation spectrum, and we calculate density 
profiles and "breathing" modes in the harmonic trap. 

We use the asymptotic Bethe Ansatz (BA) [l8| which 
is justified at sufficiently small densities, when only two- 
particle collisions are important [l9| . The underlying 
idea goes back to the earlier days of high energy physics 
and was known as S-matrix theory [20J in the 1950s. 
The BA method considers the scattering matrix be- 
tween asymptotic states as an alternative to a Hamil- 
tonian/Lagrangean description. The scattering of ultra- 
cold atoms in ID gases close to resonance can naturally 
be described by the scattering phase shift, whereas the 
formulation of a microscopic quantum Hamiltonian is dif- 
ficult. It can be shown that the scattering matrix close 
to resonance corresponds to a highly singular, although 



local, two-body interaction in the spirit of Refs. [161. |21|. 
To avoid difficulties related to the determination of the 
operators and states in this case, we use an approach 
based entirely on the scattering phase shift. 

Let us start by briefly reviewing the 3D scattering 
properties in a p-wave channel. At low energies, the 
phase shift S p (k) can be expanded as 



k 3 cot 6 p (k) = -1/wi - ai k 2 + 0(/c 4 ), 



(1) 



where wi is the scattering volume, ot\ is the effective 
range, and k is the relative momentum. For a\W\ < 0,the 
scattering length obtained from S p (k) has a very sharp 
shape resonance at the wave vector [22l. |23| 

k r = l/y/—atiwi. (2) 

Such resonance is absent for s-wave scattering, but for the 



2 



p-wave channel it exists due to the presence of the effec- 
tive range parameter. The higher order terms in Eq. ([T]) 
do not significantly affect the shape resonance, since they 
are suppressed by powers of the small parameter kR <C I , 
where R is the characteristic radius of the 3D potential. 
For fermions, the typical momenta of scattering parti- 
cles are of the order of the Fermi momentum kp = irn. 
Therefore, the condition kR <C I necessary for neglecting 
three-particle collisions 1191 implies the low-density limit 
nR < 1. It is known [13 HI HI that a x > 1/R > 
does not change significantly at the p-wave Feshbach res- 
onance, while wi can be tuned to very large absolute 
values compared to its characteristic values of the order 
of \wi \ ~ R 3 away from the resonance. 

Under transverse harmonic confinement with fre- 
quency cj_l, only the lowest transverse mode is occupied 
if the momenta of scattering fermions satisfies 



ka± <C 1, 



(3) 



where <zj_ = \J Ti / (rmu j_) , and m is the atomic mass. Un- 
der such conditions, the ID scattering amplitude in an 
odd channel is given by [Toj ] 



j?odd 

Ip — 



-ik 



l/l p + ik + k 2 ^ p ' 



where 24 



l D — 3a 



^ - 3V2£(-l/2) 
wi 



aia\ 



(4) 



> 0, (5) 



and 3V^C(-l/2) « -0.88. The notation adopted is that 
of Ref. (Tg. Using estimates [13,111 of cti for 6 Li and 40 K 
atoms at resonances with Bo s» 215G and 198. 6G, and 
with transverse frequencies ui± = 2ir x 200kHz and 2-7T x 
30kHz fl3j| . we obtain £ p /aj_ s» 50 and ps 13, respectively. 
Thus, under typical experimental conditions needed to 
achieve the ID regime £ p 3> a± and hence all three terms 
are significant in the denominator of Eq. (j4]) . 

The many-body fermionic wave function 
ip(zi, . . . , zm) is anti-symmetric, tp(. . . , Zi, . . . , Zj, . . .) = 
—ip(...,Zj,..., Zj, . . .), and discontinuous when two coor- 
dinates coincide [16[. We define its symmetrized version 

by ij} + (zi,..., zm) = rii<j Si g n ( z i ~ z j)^( z i, • • • j z m), 
which is continuous. Then Eq. (UJ) implies the following 
boundary condition: 

lim ( r + a 2 -e p 9 2 2 ) i>+(z 1 ,...,z M )=0. (6) 

z=Zj — Zi— ^0+ \ bp J 

Solving the two-body problem as i/j+(zi, Z2) oc e lA l Zl_Z2 l, 

we obtain two roots A± = i ± — 4£ p /^ /(2£ p ). 

For l p > 0, ImA + > 0, which corresponds to a bound 
state. The lowest energy state satisfying the boundary 
condition (|6|) can then be constructed as tp(z x , . . . , zm) oc 
n 4< j Sign(zi - Zj) ll^j exp(iX + \zi - Zj\). As in the at- 
tractive LL model, its energy does not have a proper 



thermodynamic limit, we will not consider the case where 
l p > 0. For l p < we construct an exact wavefunction 
tj)+(zi, . . . , zm) as a combination of plane waves, using 
the BA method in a similar manner to the LL model. In 
our case, such construction leads to the following periodic 
boundary conditions on a circle of length L 



i\ < L 



M 

n 



£p(Aj — ^fc) 2 — j^j + i(Aj — Afc) 



1 £p(<\? ~ A fc ) 5 



i(Xj - Afc) ' 



(7) 



and the total energy is given in terms of quasimomenta 
Xi as E = h 2 / (2m) X] Af . We prove that all solutions of 
Eq. ([7]) are real by writing the k-th term in the product 



as 



-A-) 



and Z p < 0, we then have 



Since ImA± < for £ p > 
< 1(> 1) for 



' (A-A+)(A-A_) 

(a-a;xa-a*) 

ImA < 0(> 0). After that, the proof simply follows the 
steps for the LL model described on p. 11 of Ref. prij . 

To obtain a thermodynamic limit, we take a logarithm 
of Eq. (0, which is written as LXj + Y^k=i ^(Aj — Afc) = 
27rnj, where rij are integer quantum numbers for odd 
M. The phase shift 0(X) is a monotonic antisymmetric 
function defined by 



0(A) = 2Arg(U - £ P A 2 + l/|y ), 



(8) 



and belongs to the interval (— 2ir, 2ir), unlike the LL 
phase shift, which belongs to the interval (— 7r, 7r). We 
then directly follow Ref. 27) and show that real solutions 
of the B A equations exist for any choice of quantum num- 
bers rij . Their values for the ground state can be fixed by 
comparison with the LL model 0, US H3| , and are given 
by rij = j — (M + l)/2. Introducing a positive function 



K(\,(i) = 9'(X,fx) 



2ig [1 + M£ p (a - m) 2 ] 
[i-\i p \^x-^] 2 +B(x-^y 



we pass to the thermodynamic limit, and write an equa- 
tion for the ground state quasimomenta distribution in 
the usual way 27rp(^) — J_ g K{y, n)p(/i)dn = 1, where 
±q is the highest (lowest) filled quasimomentum and the 
normalization is given by n = M/L = J"* p(y)dv . Apart 
from new definitions of 9(X) and K(X,n), the structure 
of the theory is similar to the LL model, and we can 
study the ground state energy, excitation spectra and fi- 
nite temperature properties using standard methods [26| . 

We choose two dimensionless parameters that deter- 
mine the ground state properties in the stable region 



7i 



IpTl 



> 0, 72 



> 0. 



(9) 



In Fig. [T]we show the dimensionless ground state energy 
functional 6(71,72) obtained by numerically solving the 
equations for the ground state. The dimensionless form 
is given by the expression 



£/L = e( 7l ,7 2 )(?m) 2 /(2m), 



(10) 
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FIG. 1: (Color online) Energy functional 6(71,72) as a func- 
tion of dimensionless parameters 71 and 72, see Eqs. (f9)l- (|10[) . 
Arrows indicate the regimes of weakly interacting fermions, 
Cheon-Shigehara (CS) model [16:], fermionic Tonks- Girardeau 
(fTG) gas [17j. and super- fTG regime (see Eq. (fTTj) ') existing 
due to the finite effective range £ p . 

and reduces to the LL functional e(7i) for 72 ^ 1, since 
the CS model [l6| obtained in this limit is dual to the 
LL model. The function 6(71,72) equals if 71 = or 
72 = 0. Expansion by methods of Ref. [28j at 71,72 3> 1 
yields e(7i,7 2 ) « e(7i) - 327r 4 /(157i7 2 ). 

We use the function 6(71,72) obtained numerically to 
evaluate density profiles in a harmonic trap within the 
local density approximation 29]. We also use the func- 
tion 6(71,72) to find the "breathing" mode frequency u> 
by solving the hydrodynamic equations 0, [3fj| . The re- 
sults are shown in Fig. [2] and depend on two dimension- 
less parameters, 71 (0) and 72(0), in the center of the 
cloud. The presence of the effective range strongly af- 
fects the shape of the profile compared to the CS model 
if 72(0) is small and 71 (0) is not too large. This effect 
can be understood by using the expansion of e(7i, 72) for 
7i j 72 1 • The leading term in the Taylor expansion 
gives 6(71,72) oc 7172 oc 1/n 2 and hence via Eq. fj 10|) the 
energy per particle for this term does not depend on den- 
sity, i.e. the gas has a divergent compressibility. Higher 
order terms in the expansion lead to a finite but large 
compressibility, which decreases with increasing 71 and a 
constant ratio 72/71- Thus, the density profile exhibits a 
strong peak near the center where 71 is smallest. 

Since fermions are more strongly correlated for 72 <C 1 
and not too large 71 compared to the fTG regime, we sug- 
gest calling such a regime a super-fTG gas, analogously 
to the super-TG gas of bosons [3l| . It is realized if 



Max( 7 2,^J)<l. 



(11) 



For \wi\ -C a.j_i the second condition corresponds to 
k r < 2nn = 2kp. Thus, the super-fTG regime is real- 
ized if the largest relative momentum of non-interacting 
fermions approaches the shape resonance wave vector 
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FIG. 2: (Color online) Density profiles under harmonic con- 
finement, measured in units of Thomas- Fermi radius Ro and 
central density no for the same cloud in the absence of inter- 
actions. Different curves correspond to {71(0), 72(0)} in the 
center of the cloud {7.5, 0.023}(blue, dotted), {97, 0.058}(red, 
dashed), and free fermions (black, thick). For a cloud of 
N = 100 6 Li atoms with uj ± = 2n x 200 kHz and uj z = 
2tv x 200 Hz, we get Ro = and these curves correspond 

to -mi/tti = 0.05,0.01,0, respectively. The inset shows the 
ratio of squares of the "breathing" and dipole mode frequen- 
cies (ui/u> z ) 2 for the same trap parameters. If a significant 
part of the cloud is in the super-fTG regime (see text), this 
ratio drops below the critical value 3 obtained in the fTG 
limit. A critical value — u)i/aj_ ~ 0.008 corresponds to a de- 
tuning from the resonance AB w 50mG. 



k r . In a similar manner to the super-TG gas of bosons, 
the super-fTG regime can be experimentally identified 
by measuring the ratio of the squares of the "breath- 
ing" and dipole mode frequencies. In the CS model 
such a ratio is always larger than 3, similar to the LL 
model [30], while the inset of Fig. [2] shows the regime 
where it is smaller than 3. A sharp decrease in this ratio 
for the super-fTG regime can be easily understood from 
the "sum rule" approach of Ref. [3(| , since the cloud den- 
sity is much more centered in the super-fTG regime than 
in the fTG regime. We can analytically estimate the 
value of —w\/a\ at which the center of the cloud enters 
the super-fTG regime, and the drop in (cu/u! z ) 2 occurs. 
For that we use Eq. Ijlip with the free fermion density in 



the center and obtain 



super-fTG 



24AT{ p / Qj 



which 



gives 0.008 for the parameters seen in Fig. [2] 

The excitation spectrum e(k) in a uniform cloud is also 
significantly modified in the super-fTG regime compared 
to predictions of the CS model. In Fig. [3] we illustrate 
several qualitative features, which appear due to the fi- 
nite effective range of interactions. Firstly, the system 
has a regime where the energy of the particle-like excita- 
tion is smaller than the energy of the hole-like excitation. 
Since the energy of the particle-like excitation should ap- 
proach fc 2 /(2m) at high momenta, there should also be 
an energy crossing. This crossing will manifest itself as 
a kink in the k— dependence of the lowest energy of the 
density wave excitations. Secondly, the spectrum of hole- 
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e(k) 




0.5 1.0 1.5 2.0 

FIG. 3: (Color online) Excitation spectrum for 71 = 0.5, 72 = 

0.005, where momentum is measured in units of ItF = 7rn, 
and energy is normalized to give a unit velocity v = e'(k)\k=o- 
Due to interactions, the particle-like mode (red, dashed) has 
energy lower than the hole-like mode (blue, solid) for suffi- 
ciently small momenta. Velocities of two modes at k = 
coincide. The inset shows energies of hole-like modes for 
71 = 10.8, 72 = 0.022 (blue, solid) and 71 = 27.3, 72 = 0.055 
(red, dashed), normalized to their respective velocities. 

like excitations can have a "roton" minimum (or even an 
additional maximum, see inset in Fig. [3]) at k « fcp. This 
minimum can be understood as a tendency of the system 
towards pairing when the parameters 71 , 72 approach the 
boundary of the stable region. Indeed, the energy of the 
hole-like excitation vanishes for k = 2kp = 2nn. Since 
particle density is twice the density of pairs, the vicin- 
ity of the paired region manifests itself as a soft mode 
at k ~ 2itn/2 = kp. Dynamic response functions of the 
system will have power-law divergences at the particle- 
and hole-like modes, which can be calculated using the 
methods of Ref. 32j. Note however, that the existence 
of the roton minimum and the inversion of the particlc- 
and hole-like spectra lead to modifications of the phe- 
nomenology of Ref. (3^1 ■ 
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